Has cosmological dark matter been observed? 
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0\ . There are many indications that ordinary matter represents only a 

0\ ■ 

tiny fraction of the matter content of the Universe, with the remainder 
assumed to consist of some different type of matter, which, for various 
■^j- ■ reasons must be nonluminous (dark matter). Among these indications 

are the inflationary scenarios which predicts that the average energy den- 
sity of the Universe coincides with the so called critical value (for which 
the expansion never stops but the rate of expansion approaches zero at 
very late times). At the same time it is known (from the predictions of 
, Big Bang nucleosynthesis on the abundances of the light elements, other 

than Helium) that the baryonic energy density (ordinary matter) must 
represent (1.5 ± 0.5)/i~ 2 % (where h is the Hubble constant in units of 
100 km s -1 Mpc -1 ) of this critical value fl],§]. We present here evidence 
supporting the model in which the rest of the energy density corresponds 
to a scalar field, which can be observed, however indirectly, in the os- 
cillation of the effective gravitational constant, and manifests itself in 
the known periodicity of the number distribution of galaxies [|3|,^]]. We 
analyze this model numerically and show that, the requirement that the 
model satisfy the bounds of light element abundances in the Universe, as 
predicted by Big Bang nucleosynthesis, yields a specific value for the red- 
shift-galactic-count oscillation amplitude compatible with that required 
to explain the oscillations described above [^|J6|, and, furthermore, yields 
a value for the age of the Universe compatible with standard bounds 
[B|. The fact that the model has successfully passed these tests lends 



support to it and to the conclusion that « 98% of the energy density 
of the Universe is stored in this scalar field. This leads to the aston- 
ishing conclusion that if the observations of Broadhurst et. al. []3|,|| are 
not a statistical fluke, the cosmological component of dark matter might 
already have been observed. 
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It has long been suspected that ordinary matter (photons, ordinary massless neutrinos 
and atoms) represents only a tiny fraction of the matter content of the Universe. The 
photons and neutrinos (collectively called radiation, because they are massless particles) 
represent at the present time a negligible fraction of the energy density corresponding to 
the atoms. The latter is sometimes called baryonic matter because most of the energy 
corresponds to the rest mass of protons and neutrons (baryons). The evidence that there 
is much more matter in the Universe comes from various sources, and it manifests itself 
through its gravitational influence. One such item of evidence is found in the dynamics of 
the motion of stars in galaxies, where the use of the standard gravitational laws to explain 
rotation curves results in the conclusion that the amount of mass contained in the galaxy 
is about ten times larger than the sum of the masses of the stars, gases and other known 
luminous constituents of the galaxy 0. In the behavior of clusters of galaxies there are 
further indications that there exists much more mass than we observe in form of stars and 
gases || . Finally, the most dramatic indication of the "insignificance" of ordinary matter 
relative to the total matter content of the Universe comes from the so-called inflationary 
scenarios P JT0|JTT[] . This indication consists in a concrete prediction that the sum of all 
contributions to the average energy density of the Universe exactly equals the so-called 
critical value. These inflationary scenarios that are invoked to solve serious defects of 
the standard cosmological model indicate that, since the observed ordinary matter is 
of the order of 1 % of the critical value, there must be a large amount of dark matter. 
Moreover, this dark matter must be mostly exotic matter, since cannot consist of radiation 
or baryonic matter. This conclusion results from the dependence on the baryonic energy 
density, of the predicted values of the relative primordial abundances of light elements 
(Be, Li, D, etc) in Big Bang nucleosynthesis. This restricts the value of the baryonic 
energy density f2b ar to lie in the range (1.5 ± 0.5)h~ 2 % of the critical value |I],|2|], where 
h is the Hubble constant in units of 100 km s Mpc . 

Cosmologists and particle physicists have long been puzzeled about what this exotic 
matter might be |12[ . The main models can be divided in two categories. The first one is 
called hot dark matter which consists of light neutrinos or a similar species, i.e. massive 
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particles whose number density at freze out, is determined when they still may be con- 
sidered as relativistic particles. The cold dark matter model consists of all the remaining 
weakly interacting massive particles like axions, neutralinos, superheavy monopoles, pri- 
mordial black holes, etc. These particles were already non relativistic when their number 
density reached a state of equilibrium in which annihilations freezes out. No independent 
evidence for the existence of these new types of matter has so far been found. 

An apparently completely independent problem in cosmology is posed by the recent 
observations in deep pencil beam surveys @,(|] which show that the number distribution 
of galaxies exhibits a remarkable periodicity. This is a shocking development, since if 
taken at face value it would imply that we live in the middle of a pattern consisting of 
concentric two-spheres that mark the maxima of the galaxy number density. This, of 
course, would lead to catastrophic consequences for our concepts of cosmology. While it 
is true that such periodicity has been observed only in the few directions that have been 
explored so far, it would be a remarkable coincidence if it turns out that it is absent in 
other directions, and we have just happened to have chosen to explore the only directions 
in which this phenomenon is observed. It seems, therefore, reasonable to assume that 
the periodicity will be also present in deep pencil beam surveys in other directions, thus 
forcing us toward the concentric spheres scenario. The seriousness of the situation is 
such that this type of scenario has indeed been proposed in a model where the formation 
of these concentric shells is a result of a "spontaneous breakdown of the cosmological 
principle" via a mechanism that results in the appearance of patches filled with a pattern 
of concentric spheres, with these patches filling the Universe It would be difficult 

to explain how we we managed to be living in the center of such a patch (more precisely 
inside the innermost sphere of one such patch). 

The only known escape from this type of scenario is to assume that the spatial peri- 
odicity is only an "illusion" and is the result of a true temporal periodicity which affects 
our observations of distant points in the Universe, and which is mistakenly interpreted as 
spatial periodicity [J5J. The models that have been put forward that allow this temporal 
periodicity involve the oscillation of an effective coupling constant due to a contribution 
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from the expectation value of some scalar field that actually oscillates coherently in cos- 
mic time at the bottom of its effective potential. It is worth to mention that a completely 
different scenario consisting of oscillating peculiar velocities of galaxies have also been 
proposed for explaining the observed periodic redshifts W% . 

The specific oscillating-coupling-constant models which have been proposed involve 
the oscillation of the effective electric charge, electron mass, galactic luminosity or the 
gravitational constant PJ6| JT5|JT6|fi7|] . Of these, the first two have been shown to conflict 
with bounds arising from tests of the Equivalence Principle |Tj|]. The model of oscillating 
galactic luminosity is a vaguely specified scenario which nevertheless seems to require at 
least two new hypotheses: a new type of star cooling mechanism, and also (as the other 
alternatives do) an oscillating cosmological scalar field which turns that mechanism on 
and off periodically in cosmic time. Needless to say, once the scenario is implemented in 
a specific model, unforeseen new bounds might also have to be overcome. 

In this light the Oscillating G Model seems to be the most attractive alternative. We 



also should point out that in the fossil record of marine bivalve shells [Tj| , there seems to be 
further evidence in support of an effective gravitational constant that oscillates with time. 
Moreover, cosmological consequences of high-frequency oscillations of Newton's constant 



have been also studied [19|. In this work, we consider the model, initially proposed in 

||, and further studied in |^,|r^,|T^,|T^] , of a massive scalar field non-minimally coupled to 
gravity leading to an effective gravitational constant which oscillates in cosmic time and 
in this way produces the illusory spatial periodicity observed in the number distribution 
of galaxies. 

The central feature of this model is a cosmological massive scalar field <p that is non- 
minimally coupled to the curvature of spacetime. The corresponding Lagrangian can be 
written as [HI 
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where Go is the Newtonian gravitational constant, £ is the non- minimal coupling constant, 
m is the mass associated with the scalar field, and £ ma t represents a matter Lagrangian. 
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The non-minimal coupling of the scalar field <p to curvature results in an effective gravita- 
tional "constant" G e ff = GV(1 + 167rGo£0 2 ) which explicitly depends on the cosmic time 
due to the contribution of the expectation value of the scalar field. 

We studied the resulting cosmological model corresponding to an isotropic and ho- 
mogeneous spacetime; i.e., the one described by the Friedman- Robertson- Walker (FRW) 
metric. It was also assumed that the scalar field, as well as the matter fields, possess the 
same symmetries as the spacetime. 

We studied the model in detail using numerical methods, and explored its compatibility 
with standard cosmological tests. The model is known to require a "fortunate phase" 
in order to satisfy the bounds imposed on the value of G/G by the Viking radar echo 
experiments p0[|5[, and by the limits on the Brans-Dicke parameter |J. It has also been 
argued || that it requires fine tuning of the parameters and data in order to satisfy the 
bounds of Big Bang nucleosynthesis and therefore that the model should not been taken 
seriously. 

We will argue that, if seen in the proper context, this fine-tunning can be taken 
instead, to be a concrete prediction of the values of parameters and more specifically as 
relationships among them, a prediction that can be used to rule out or support the model. 

Our general philosophy concerning this fine-tuning should be understood in the fol- 
lowing context. Scientific models that require a very precise choice of the numerical value 
of the initial conditions in order to reproduce a given qualitative behavior of the observa- 
tional data are models that would be considered unnatural, and the choice of the specific 
initial data is justifiably described as "fine-tuning" . However, models that require a very 
precise choice of the numerical value of the initial conditions in order to reproduce a 
specific numerical observational data cannot be considered as unnatural, especially if for 
every conceivable value of the observational data (at least in some range), there is a corre- 
sponding value of the initial data. In this type of models, the particular "preferred" value 
of the initial data is just the result of a 1 to 1 correspondence between initial conditions 
and final outcome. The exactness of the predictions so obtained will, of course, depend on 
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the exactness of the observational data that determine the remaining parameters entering 
in the model. 

The main point is that in order to reproduce a specific value for the red-shift-galactic- 
count oscillation amplitude Aq and the corresponding periodicity of 128 Mpc h~ x , a very 
specific value for the density of the baryonic matter in the Universe emerges. In a previous 
work fL7l , for example, we took the quoted observational value of Ao = 0.5 [|J and h = 1 
for the Hubble constant. The resulting value, £\ar ~ 0.021 for which it was possible to 
recover the primordial 4 He abundance from nucleosynthesis and which was also consistent 
with the inflationary value Q = 1 lay in a very narrow range that we had assumed to be 
0.016 < fibar < 0.026 [[Hp and thus we took it to be a strong piece of evidence in favor of 
the model. Alternatively, if we start from a specific value of the baryonic energy density 
(for example one that lies within the range allowed by nucleosynthesis) and the frequency 
of the periodicity, the result is a specific value for the amplitude Aq. Moreover, these 
values, together with the value of the Hubble parameter, determine a specific value for 
the age of the Universe, and so it is a highly nontrivial question whether these results are 
or are not compatible with observations. Our main result is to answer the above in the 
affirmative, thus implying that the problems of the nature of cosmological dark matter 
and of the periodicity in the galaxy number distribution may be solved simultaneously 
within the framework of an oscillating G model with a massive scalar field. 

The field equations following from the Lagrangian ([]]) are similar to Einstein's equa- 
tions with the Newtonian constant replaced by the effective gravitational constant G e Q, 
and an energy-momentum tensor which contains contributions from the non-minimal cou- 
pling to curvature, the scalar field, and the matter. The last one consists of two non- 
interacting perfect fluids: TZt = ?C + TJf = £*=i,2 [(pi + e i )W i U v + p l g fiu ] . The first 
one corresponds to pure baryonic matter (i = 1), and the second one represents a pure 
radiation field (i — 2). Moreover, the scalar field must satisfy the massive Klein-Gordon 
equation with an extra term due to the non-minimal coupling to curvature. 

Our analysis consists in evolving the scale factor, the scalar field and the ordinary 
matter densities backwards and forwards in cosmic time using the field equations, and 
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starting from the model parameters £ and m, and data corresponding to today's values 
of H , ^bar,^rad, 0o and 0o- With these initial conditions it is then possible to integrate 
the field equations numerically. 

We will work within the standard inflationary scenario, so Q = f2bar + ^rad + = 1, 
where f2 ra d is the cosmic background radiation (CBR) and its value today is determined 
from the CBR temperature of 2.735 K resulting in a contribution which is several orders 
of magnitude smaller than Q^ aT which is itself known from recent analysis to be in the 
range 0.01/i~ 2 — 0.02/i -2 0,0] ■ Therefore fixing Q^ av is equivalent to fixing Q^. 

The initial condition (today) for the time derivative of the scalar field <po = was 
fixed so that it satisfies the Viking radar echo experiments [50,51. The initial condition 
4>q as well as the value of the coupling constant £ turn out to be expressible in terms 
of the values of the amplitude Aq, the oscillation frequency uj of the scalar field and the 
parameter = 1 — f2bar — ^rad- The observations yield a value of Aq of about 0.5 (in fact, 
it has been argued that Aq > (9(0.5) ||, see the discussion below), while the observed 
galatic periodicity of 128 h -1 translates to a value m ~ 10 -31 eV for the scalar-field mass 

» 

As we mentioned, one of the main problems faced by the model is that related to 
the primordial nucleosynthesis of 4 He. The latter is determined by the temperature at 
which the rate of weak interactions, which convert neutrons to protons, equals the value 
of the Hubble parameter. According to standard cosmology the value of about 0.7MeV 
for this temperature corresponds to a 4 He abundance which is the best approximation to 
observational data (see P^j for a review). 

The evolution of the scalar field backwards in cosmic time results in it going to ±oo 
depending on the initial data. This suggests that there is a very precise initial data for 
which it is possible to reach a transition point in which remains steady and close to 
zero. This "steady state" is represented by a kind of plateau during which G c g — > Gq 



17 1 - It was possible to correlate the length of this plateau with the recovery of the 



precise standard freeze-out temperature. The larger the plateau the closer the freeze-out 
temperature predicted by the oscillating model approached the value 0.7 MeV. The search 
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for that transition point is what we call "fine- tuning' 1 and it turns out that this can be 
done by adjusting only the parameters Aq and f2bar (see Fig. [[]). In principle it is possible 
to extend the plateau of <fi to the nucleosynthesis era or even to earlier eras by improving 
the "fine- tuning" of the values of fl^ ar or Aq. 

As explained above, the recovering of the success of standard Big Bang nucleosynthesis 
requires the implementation of the so-called "fine-tuning 11 of the parameters and the initial 
conditions of the model. However, as we previouly argued this is not the kind of fine- 
tuning that implies the dismissal of the model, but it is rather a procedure that becomes 
necessary in order to recover the observational data extracted from our Universe today. 
We must also stress that while the fine tuning is completely unnatural when approached, 
as we have approached to it, from the present to the past, when looked from the opposite, 
and more natural direction, the situation is quite different. In fact all that seems to be 
required is for some mechanism to drive the scalar field to an extremely low value before 
the era of nucleosynthesis. Then, as our calculations show, the field will remain at that 
value up to and beyond that era so that we will have G e g ~ Go and then the success of 
"Big Bang Nucleosynthesis" will be recovered naturally. The value of the field will later 
be amplified by the curvature coupling just before the onset of oscillatory behavior. 

The era with 0^0 ends just before the Hubble parameter approaches to the value of 
m, a situation that is followed by an amplification and then the onset of oscillations of <f). 

But the point is that this would actually represent no "fine-tuning" at all (if looked in 
the right perspective), because all that it will mean is that (given the physical constant 
precise values) the standard inflationary mechanism will ensure that the energy densities 
of the various components, i.e., scalar field, baryon matter and radiation, add up to 
Q — 1, which will then correspond then to a Universe at our time with precise values 
of the densities, expansion rate etc. In particular the precise current value of the scalar 
field amplitude and phase arise from a particular precise value of the parameters at early 
times, among them the baryon content of the Universe. 

Thus it is possible that starting from an arbitrary value of <fi near the Big Bang, a 
mechanism related to inflation would drive the scalar field to a value near zero, where 
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it will remain until just before H m m when amplification and then oscillations would 
occur. Now, since the behavior of G e g is determined by the expectation value of <j), 



oscillations of <fi induce oscillations in G c s ||17|| , and so we are led to the scenario in which 
the temporal oscillations of the effective gravitational constant manifest themselves in an 
apparent spatial periodicity of the number distribution of galaxies [|5,15,17 . 



As we have said, our initial research was carried out by taking very specific values of 
the observed quantities The present analysis is motivated by our previous results, 
and by the need to study the compatibility of the model given the uncertainties in the 
most up-dated observational data. In particular we consider the constraints arising from 
the lower bound on the age of the Universe for which we take the most conservative 
lower bound of 11.5 billion years obtained from a giant-branch fitting [fJJ. In addition, 
observations of the galaxy number distribution do not fix the amplitude, but instead 
indicate the lower bound Ao > (9(0.5). Finally, for the Hubble parameter we will use 
the range of values 0.65 < h < 0.75 which is the intersection of the ranges allowed by 
recent astronomical observations (data from Hubble Space Telescope, and studies of type I 
supernovae p^ , |23| fl|). It is to be emphasized that the largest uncertainty in the observed 
galactic periodicity of 128/i -1 is that contained in the value of h. Therefore, we will 
investigate the sensitivity of our results with respect to the values of Ao, and h. 

For each value of h in the above range we obtain a corresponding allowed range of 
Qbar, and for each value of f2bar we can obtain by means of the "fine-tuning" mechanism 
a value of Ao and a resulting value for the age of the Universe. We have performed this 
analysis with h = 0.65 and h = 0.75 the extreme values of the interval. In the first case, 
the range allowed by light-element nucleosynthesis corresponds to 0.023 < Qbar < 0.047, 
the corresponding range in Ao was found to be [0.490 , 0.499] and the age of the Universe 
lies in the range [11.77, 12.06] Gy which satisfies the 11.5 Gy bound. In the second case 
the range allowed by light-element nucleosynthesis corresponds to 0.023 < Qbar < 0.046, 
the corresponding range in Ao is [0.4950 , 0.5015] and the age of the Universe lies in the 
range [10.32 , 10.57] Gy which fails to satisfy the 11.5 Gy bound. 

Similar results can be obtained for other values of h contained in the range 0.65 < 
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h < 0.75, from which we find that the upper value of h that satisfies the age bound is 
h ~ 0.68 

These results show that there is a range of the parameters of the model that satisfies 
all known cosmological constraints while at the same time explaining both the nature of 
the cosmological dark matter and the observed periodicity in the galaxy distribution with 
red-shift. Figure || shows the red-shift periodicity of the Hubble parameter as calculated 
by the oscillating G model. 

Figure ^| shows the bounds imposed on J\ ar with h = 0.65 and h = 0.75, and on the 
age of the Universe from a giant-branch fitting 0. A triangle indicates a configuration 
calculated with our model with Q ~ 0.0336 and A = 0.495. Figure f| shows the freeze-out 
temperature predicted by the same configuration (~ 0.703 MeV). The age of the Universe 
extracted from data of Fig. |5| corresponds to 0.792//^ 1 which translates into 11.917 Gy 
with h = 0.65. 

The cross in Fig. ^| represents a configuration with h = 0.75 and Q ~ 0.0210 and 
A = 0.5. Although the results are compatible with almost all bounds, the resulting age 
of the Universe 10.505 Gy turns to be too small and thus this configuration must be ruled 
out. 

To conclude, the oscillating G model is certainly the most attractive model for explain- 
ing the observed periodicity in the galactic distribution, and it should also be considered 
as a missing mass model with the scalar field playing the role of cosmological dark matter 
which is, however, indirectly observable in the oscillation of the galactic distribution. The 
requirement that the scalar field go through a plateau phase during which G e g ~ Go then 
yields a relationship between the value of the baryonic energy density and the value of the 
redshift-galaxy-count oscillation amplitude. The fact that these can be accommodated 
within the allowed ranges inferred from observations, together with the fact that they 
lead to an acceptable value for the age of the Universe is taken as support for the model, 
and thus for the idea that the rest of the energy density of the Universe is stored in the 
oscillating scalar field. This leads to the astonishing conclusion that, if the observations 
of Szalay et. al. are not a statistical fluke, the cosmological component of the dark matter 
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might have been observed. 

One remaining problem to be faced is that of the phases which according to the latest 
observations are not the same in all the explored directions. One possible way out of this 
difficulty has been suggested in |J where it is argued that a modification of the form of 
the scalar potential helps to ameliorate the problem. A careful analysis of such a scenario 
should probably be done within the context of the analysis of local inhomogeneities. 
Moreover the model should of course be extensively tested, and in particular further 
and more precise observations will be needed before one could claim with any degree of 
certainty that indeed cosmological dark matter has been observed. 
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FIGURES 

FIG. 1. The fine-tuned scalar field amplitude as a function of ln[a/ao] for a flat Universe 
(f2 = 1) with Obar ~ 0.033 689 and Ao = 0.495 at the onset of oscillations. At present time 
(a = 0) the initial amplitude is 0o ~ 3.26 x 10~ 3 and £ ~ 6.2859. Computations were stopped 
at a ~ 1.5 . 

FIG. 2. Periodic distribution of the Hubble parameter with respect the red-shift. The disti- 
bution extends from z ~ —0.5 to z ~ 0.5. 

FIG. 3. Bounds imposed by observations on the age of the Universe (horizontal dashed 
line), on the values of h and therefore on the resulting allowed ranges of the baryonic component 
of matter (vertical lines) obtained from 0.01 < Slbar h 2 < 0.02. The triangle depicts the 
configuration corresponding to Figs. [I], 0, [|, ^. The cross represents a configuration with h = 0.75 
(see text). The square indicates a configuration barely allowed by the age-of-the-Universe bound. 
This corresponds to the values h = 0.68, S\ar ~ 0.0222 and *4o = 0.4995 . 

FIG. 4. Expansion (solid line) and neutron-to-proton-weak-interaction transition (dashed 
line) rates in terms of the blackbody temperature. The asterisk depicts the freezeout temperature 
~ 0.7 MeV at which nucleosynthesis takes place as predicted by the standard cosmological 
models. The crossing point of the curves indicates the corresponding freezeout temperature 
~ 0.7 MeV for the oscillating model of previous figures. 

FIG. 5. The scale factor of the oscillating model of the previous figures in units of its value 
today as a function of cosmic time (in units of H^ 1 ). The present time to has been taken to be 
t = lH \ 



15 



